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Abstract 



We study gauge and Yukawa coupling unification in the simplest extension of 
the Minimal Supersymmetric Standard Model (MSSM) which incorporates R— Parity 
violation through a bilinear superpotential term. Contrary to what happens in the 

H , 

MSSM, we show that bottom-tau unification at the scale Mqut where the gauge 
couplings unify can be achieved for any value of tan/3 by choosing appropriately 
the sneutrino vacuum expectation value. In addition, we show that bottom-tau-top 
unification occurs in a slightly wider tan (5 range than in the MSSM. 



The Standard Model (SM) of particle physics is very successful in describing the 
interactions of the elementary particles, except possibly neutrinos. Although it is re- 
garded as a good low-energy effective theory, the SM has many unanswered questions 
and theoretical problems. Its gauge symmetry group is the direct product of three groups 
££7(3) x SU(2) x £7(1) and the corresponding gauge couplings are unrelated. It does not 
explain the three family structure of quarks and leptons, and their masses are fixed by 
arbitrary Yukawa couplings, with neutrinos being prevented from having mass. The Higgs 
sector, responsible for the electroweak symmetry breaking and for the fermion masses, has 
not been verified experimentally and the Higgs mass is unstable under radiative correc- 
tions. As a result, say, the hierarchy between the electroweak scale and the Planck scale 
is not understood. 

In supersymmetry (SUSY) |J the Higgs mass is stabilized under radiative corrections 
because the loops containing standard particles is partially cancelled by the contributions 
from loops containing supersymmetric particles. If we add to the Minimal Supersymmetric 
Standard Model (MSSM) § the notion of Grand Unified Theory (GUT), then we find 
that the three gauge couplings unify at a certain scale Mqut @- Indeed, measurements 
of the gauge couplings at the CERN e + e~ collider LEP and neutral current data [|J are 
in good agreement with the MSSM-GUT with the SUSY scale M SU sy ~ 1 TeV §. In 
addition, the unification scale in SUSY-GUT is high enough to predict a proton decay 
rate slower than present experimental limits, as opposed to the non-SUSY GUTs, where 
the proton decays too fast. 

Besides achieving gauge coupling unification |J, GUT theories also reduce the num- 
ber of free parameters in the Yukawa sector. For example, in £77(5) models, the bottom 
quark and the tau lepton Yukawa couplings are equal at the unification scale, and the 
predicted ratio mb/m T at the weak scale agrees with experiments. Furthermore, a relation 
between the top quark mass and tan /3, the ratio between the vacuum expectation values 
of the two Higgs doublets is predicted. Two solutions are possible, characterized by low 
and high values of tan/3 [J7|. In models with larger groups, such as £0(10) and E 6 , both 
the top and bottom Yukawa couplings are unified with the tau Yukawa at the unification 
scale M. In this case, only the large tan/? solution survives. 



Recent global fits of low energy data to minimal supersymmetry QTOj show that it 



is hard to reconcile these constraints with the large tan /3 solution. Specially important 
are the measurements of the B(b — > 57) decay rate and the bound on the lightest Higgs 
mass. In addition, the low tan f3 solution with fi < is also disfavoured. In this letter, we 
show that the minimal extension of the MSSM-GUT [llj i n which R-Parity violation is 
introduced via a bilinear term in the MSSM superpotential |12|, [13[ , allows b — r Yukawa 
unification for any value of tan/3 = v u /vd and satisfying perturbativity of the couplings. 
We also analize the t — b — r Yukawa unification and find that it is easier to achieve than 
in the MSSM, occurring in a slightly wider high tan (3 region. 
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For simplicity, we consider only the third generation of quarks and leptons. In this 
way, the superpotential is given by 



W = h t Q 3 U 3 H u + h b Q 3 D 3 H d + h T L 3 R 3 H d + fiH u H d + e 3 L 3 H u 



where the first four terms correspond to the MSSM and the last one is the bilinear term 
which violates R-parity. This superpotential is motivated by models of spontaneous 
breaking of R-Parity [H|. Here, R-Parity and lepton number are violated explicitly by 
the e 3 term. 

It is clear from eq. ([I]) that the scalar potential contains terms which induce a non- 
zero vacuum expectation value (VEV) of the tau sneutrino (v T ) = v 3 / y/2. It contributes to 
the W mass according to m^y = \g 2 (v d + v\ + v 3 ), where v d j a/2 and v u j a/2 are the VEVs 
of the two Higgs doublets H d and H u respectively. The R-Parity violating parameters 
e 3 and v 3 violate tau-lepton number, inducing a non-zero v T mass m UT oc (/j i v 3 + e 3 v d ) 2 , 
which arises due to mixing between the weak eigenstate v T and the neutralinos. The 
latest v T mass limit from ALEPH is iii Vt ^ 16 MeV. The v e and remain massless in 
first approximation. They acquire typically smaller masses from supersymmetric loops. 
As already mentioned, in what follows we consider only the third generation of quarks 
and leptons. 

It is important to note that the e-term in eq. ([I]) is a physical parameter and cannot 
be eliminated by a redefinition of the superfields H d and L 3 fl5f . The reason is that, 
after the rotation, bilinear terms which induce a tau sneutrino VEV are re-introduced 
in the soft scalar sector |16|]. Moreover, in contrast to many prejudices |T7|], we wish to 
stress that the R-Parity violating parameters v 3 and £3 need not be small. In models with 
universality of soft supersymmetry breaking mass parameters |16| m Vr is naturally small 
because it arises from a seesaw mechanism in which the the effective mixing arises only 



radiatively, and can easily lie in the eV range |11 



R-Parity violation also implies that the charginos mix with the tau lepton, through 
a mass matrix is given by 



M 



c 



M -4 



V5 gVd 



-£3 





^h T v d 



(2) 



with h T being the tau Yukawa coupling. Imposing that one of the eigenvalues reproduces 
the observed tau mass m T , the tau Yukawa coupling can be solved exactly as |L3| 



2m 2 

v d 



l + 5 2 



(3) 



where the 5; 



1,2, depend on m T , on the SUSY parameters M, /x, tan/3 and on the 



R-parity violating parameters €3 and t^. They can be found in ref. [13[ and can easily be 
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shown to vanish in the MSSM limit e 3 — > and V3 — > 0. On the other hand, the bottom 
and top Yukawa couplings are related to the bottom and top masses according to 

v v 
m t = h t — y= sin /3 sin # , m b = h b — y= cos/? sin 9 (4) 

where we use spherical coordinates for the VEVs, defining v = 2m w /g, tan/5 = v u /vd, 
and cos 9 = v 3 /v. 

We now turn to the study of the renormalization group evolution of the various 
relevant parameters of the model such as the gauge and Yukawa couplings, the SM quartic 
Higgs coupling and the third generation fermion masses. In our approach we divide the 
evolution into three ranges: (i) from M S usy to M GUT , where we use the two-loop RGEs 
of our model, (ii) from m t to Msusy, where we use the two-loop SM RGEs including the 
quartic Higgs coupling and (iii) from Mz to m t we use running fermion masses and gauge 
couplings. 

Using a top-bottom approach we randomly vary the unification scale Mqut and 
the unified coupling olqut looking for solutions compatible with the low energy data |1| 
a-^{m z ) = 128.896 ±0.090, sin 2 9 w {m z ) = 0.2322 ± 0.0010, and a s {m z ) = 0.118 ±0.003. 
We use the approximation of a common decoupling scale Msusy ~ 1 TeV for all the 
supersymmetric particles. The solutions we find are concentrated in a region of the 
Mqut — &gut plane. For the simpler case where the SUSY scale coincides with the top 
mass, Msusy = mt, this region is centered at the point Mqut ~ 2.3 x 10 16 GeV and 
acuT^ 1 ~ 24.5, which we adopt from now on. 

Next, we study the unification of Yukawa couplings using two- loop RGEs. We take 
m w = 80.41 ±0.09 GeV, m T = 1777.0 ±0.3 MeV, and m b (m b ) = 4.1 to 4.5 GeV pf. We 
calculate the running masses m T (m t ) = r q~ 1 m T (m T ) and m b (mt) = ri b ~ l m b {m b ) 1 where 7] T 
and rj b include three-loop order QCD and one-loop order QED |J. At the scale Q = m t 
we keep as a free parameter the running top quark mass m t (m t ) and vary randomly the 
SM quartic Higgs coupling A. Using SM RGEs we evolve the gauge, Yukawa, and Higgs 
couplings from Q = m t up to Q = Msusy- The initial conditions for the SM Yukawa 
couplings are A|(m t ) = 2mf{m t )/v 2 ) with i = t, b, t and v = 246.2 GeV. 

At the scale Q = Msusy, below which all SUSY particles are decoupled (including 
the heavy Higgs bosons) we impose the following boundary conditions for the quark 
Yukawa couplings 

M m susy) = h t( M susv) sin/? sin ^ 

^b(Mg USY ) = h b (M^ USY ) cos/? sin 9 (5) 

where hi denote the Yukawa couplings of our model and Aj those of the SM. Due to its 
mixing with charginos, the boundary condition for the tau Yukawa coupling is slightly 
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m l (m l )=4.3 GeV 
<X„(M Z )=0.1 1( 




Figure 1: Top quark mass as a function of tan/? for different values of the R-Parity 
violating parameter t> 3 . Bottom quark and tau lepton Yukawa couplings are unified at 
Mqut- The horizontal lines correspond to the la experimental m t determination. Points 
with t — b — t unification lie in the diagonal band at high tan j3 values. We have taken 
Msusy = m t . 

more complicated: 



K(Ms USY ) = h T (M+ USY )cosf]sm9^l^^. ( 6 ) 
Finally, the boundary condition for the quartic Higgs coupling is given by 

KMsusy) = \{(g\M+ USY )+g'\M+ USY )}(cos2Psm 2 6 + cos 2 6) 2 (7) 

The MSSM limit is obtained setting 9 —>• tt/2 i.e. V3 = 0. 

At the scale Q = Msusy we vary randomly the SUSY parameters M, /1 and tan/3, 
as well as the R-Parity violating parameter €3. The parameter V3 = vcos9 is calculated 
from eq. (|?p. Since A (or equivalently the SM Higgs mass m 2 H = 2Xv 2 ) is varied randomly, 
in practice we also scan over 9. This way, we consider all possible initial conditions for 
the RGEs at Q = M S usy, and evolve them up to the unification scale Q = M GUT . The 
solutions that satisfy b — r unification are kept. 

In Fig. D we illustrate our point by plotting the top quark mass (we always use the 
pole mass) as a function of tan (3. For simplicity we have taken Msusy = m>t but it should 
be clear that a different M S usy choice would not change qualitatively our results. Each 
selected point in our scan satisfies bottom-tau unification (to within a 1%) hb^Mcur) = 
h T {McuT) and it is placed in one of the shaded regions according to the value of \v%\. The 
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first region with t> 3 = e% = corresponds to the MSSM and sits at the top of the plot. 
Points with |u 3 | < 1 GeV fall in the region just below. The subsequent regions labelled 
by 1 < |u 3 | < 5 GeV up to \v 3 \ > 40 GeV respectively are obtained when v 3 gets higher. 
They are narrower in tan/3. Note that points with smaller v 3 values, say 1 < |i>3| < 5 
GeV fall in the region labelled as such, as well as in all previous regions, but not in the 
subsequent ones. This overlapping with the previous regions decreases as we increase \v 3 \ 
in such a way that points with \v 3 \ > 40 GeV fall almost exclusively in the last region. 

The two horizontal lines correspond to the top quark mass within a la error. In the 
MSSM limit we can see the two solutions compatible with the experimental value of the 
top quark mass, one with tan f3 ~ 1 and the other with tan (3 ~ 55-60. It is clear from the 
figure that by selecting appropriately the value of |i> 3 | we can find b — r unification for any 
tan (3 value within the perturbative region 1 < tan (3 < 62 of the Yukawa couplings. For 
1^3 1 ^ 20 GeV one has, as in the MSSM, two disconnected solutions for b — r unification, 
one with tan (3 ~ 1, and a large tan (3 range which, for intermediate v 3 can be quite broad. 
Note that for 20 < \v 3 \ < 40 GeV only the tan/3 range from 3 to 8 or so is consistent 
with the 1 a top mass measurement, for the chosen a s and mb(mb) values. Similarly, 
the \v 3 \ range above 40 GeV would be ruled out. Note that our results do not depend 
qualitatively on the definition chosen for tan f3. For example, if we define tan f3 in the 
way which is natural in the basis where the e3-term disappears from the superpotential, 
tan (3' = v u / <Jv \ + v\ we also can find b — r unification for any tan (3' value. 

We now turn to the discussion of the uncertainties of the unification program in 
this model. The general trend follows closely that of the MSSM. The dependence of our 
results on the strong coupling constant and the bottom mass running is totally analogous 
to what happens in the MSSM. Indeed, we have studied the effect of varying a s in Fig. p] 
and found that the upper bound on tan/3, which is tan/? ^ 61 for a s = 0.118, increases 
with a s and becomes tan/3<63 (59) for a s = 0.122 (0.114). On the other hand the MSSM 
region is narrower if a s increases, specially at high tan f3 values. We have verified that the 
same trend extends to the regions with large v%. Finally, we mention that the top mass 
value for which unification is achieved for any tan (3 value within the perturbative region 
increases with a s , as in the MSSM. Turning to the dependence on nib, the behaviour is 
the opposite one. In Fig. |l]we have taken m,b(mb) = 4.3 GeV. As before the value of tan/? 
is bounded from above by tan/3 ^ 61 due to the perturbativity condition of the bottom 
quark Yukawa coupling. If we consider mbijrib) = 4.1 (4.5) GeV then the upper bound 
of this parameter is given by tan/3 < 64 (58). In addition, the MSSM region is narrower 
(wider) at high tan/3 compared with the mb(mb) = 4.3 GeV case shown in Fig. [I]. 

Finally we have studied the possibility of bottom-tau-top unification in our model. 
The diagonal line at high tan (3 values corresponds to points where t — b — r unification is 
achieved. Since the region with |t> 3 | < 5 GeV overlaps with the MSSM region, it follows 
that t — b — t unification is possible in this model for values of \vs\ up to about 5 GeV, 



5 



instead of 50 GeV or so, which holds in the case of bottom-tau unification. Within the 
MSSM, t—b—r unification is achieved in the range 56^ tan /3^59 with m t completely inside 
the la region. In this case, bilinear R-Parity violation does not enlarge the allowed tan j3 
region. However, at the 2a level our model allows t — b — r unification for 54 ^ tan (3 ^ 59. 
In addition, we have checked that the region with t — b — r unification in the MSSM case 
shrinks if a s is increased. The space left out by the MSSM is taken over by the regions 
with |t>3| < 5 GeV so that, for large a s , t — b — r unification occurs in a wider tan/3 range 
than possible in the MSSM, even in the la level. 

In conclusion, we have summarized [111 the results of the first systematic study of 
gauge and Yukawa coupling unification in a model where we introduce bilinear R-Parity 
violation. The model is the simplest alternative to the MSSM which mimics in an effective 
way many of the features of models of spontaneous breaking of R-Parity. We showed that, 
in contrast to the MSSM, where bottom-tau unification is achieved in two disconnected 
tan f3 regions, in our model b — r unification occurs for any tan j3 value, provided we choose 
appropriately the value of the tau sneutrino vacuum expectation value v 3 . In addition, 
we showed that t — b — r unification is achieved for \v^\ ^5 GeV at high values of tan/? 
in a slightly wider region than that of the MSSM. Apart from the intrinsic interest in the 
study of broken R-parity models, because of their theoretical as well as phenomenological 
importance, our results are relevant in connection with LEP bounds on the Higgs boson 
mass and recent measurements of the B(b — > S7) decay rate. Taken at face value, these 
disfavour the high tan (3 solution and also the low tan (3 solution with fi < in the MSSM 



as suggested in [10 
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